arXiv:astro-ph/9808180vl 18 Aug 1998 


Baryon Crisis and Cluster Mass Estimates 


Xiang-Ping Wu 

Institute of Astronomy, National Central University, Chung-Li Taiwan 
32054, ROC; and Beijing Astronomical Observatory, Chinese Academy 
of Sciences, Beijing 100080, China 

Abstract. I report the current status on the determinations of the 
baryon fractions and dynamical/lensing masses of galaxy clusters as well 
as the a - T relationship, making use of all the published data in litera¬ 
ture which include 304 measurements for baryon fractions, 320/21 data 
points for masses of intracluster gas/cluster galaxies, (152+228)/72 data 
points for cluster masses derived from X-ray emitting gas/optical galax¬ 
ies based on the hydrostatic equilibrium hypothesis, 46/33 measurements 
for the projected cluster masses obtained from strong/weak lensing. A 
straightforward statistical analysis yields the following features: (l)The 
cluster baryon fraction indeed shows an increasing tendency to reach a 
value of ~ 30% at large cluster radius. (2)The dynamical cluster masses 
are in fairly good agreement with the masses determined by weak lensing 
technique; (3)The cluster masses provided by strong lensing and other 
methods are essentially consistent within a factor of ~ 2. (4)The clus¬ 
ter a — T relationship is in concordance with the scenario that overall 
galaxy clusters are the well relaxed dynamical systems. These results are 
strongly suggestive of a low-mass density universe (f l m < 0.3), in which 
the fraction of baryonic matter (~ 30%) is not very small ! 


1. Introduction 

The standard Big Bang Nucleosynthesis (BBN) has set rather a tight limit on 
the baryonic matter density of the universe (Walker et al. 1991; Schramm & 
Turner 1997): 0.028 < Lib < 0.08 (I adopt Ho = 50 km s -1 Mpc -1 throughout). 
The recent measurements of various baryonic matter components in the nearby 
universe indeed yield a good agreement with the BBN prediction (Fukugita et 
al. 1997), indicative of a relatively small fraction of baryonic matter if the 
total mass density of the universe is Ll m > 0.1, i.e., a significant fraction of the 
mass in the universe is invisible. With the progress in X-ray technique over 
the past decade, a considerably large amount of hot X-ray emitting gas has 
been detected in groups and clusters of galaxies, which gives rise to a baryon 
fraction of up to fb ~ 30% (see Figure |l). If clusters are a fair sample of 
matter composition of the universe and if a flat cosmological model of flo = 1 
is acceptable, we are inevitably faced with the difficulty of the so-called baryon 
crisis (White et al. 1993). The most prevailing solution to the puzzle today 
is to work with an open universe with matter density f l m < 0.3, which is, 
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nevertheless, consistent with the result obtained by a number of recent studies 
on the abundances and evolution of galaxy clusters, in particular the mass-to- 
light ratio M/L measurements (e.g. Bahcall & Cen 1993; Bahcall et al. 1995; 
Carlberg et al. 1997; Bahcall et al. 1997). 

The key issue of either confirming or resolving the baryon crisis is closely 
connected to the question of how accurately one can determine the masses of 
galaxy clusters. Recall that the traditional cluster estimators strongly rely upon 
the assumption of hydrostatic equilibrium: both optical galaxies and intracluster 
diffuse gas trace the underlying gravitational potential of the whole cluster. An 
independent method of estimating cluster mass, which has been available only for 
decade, is to employ the gravitational lensing technique. It gives rise to cluster 
masses regardless of the cluster matter components and their dynamical states. 
A series of work has thus been made on a statistical comparison of the cluster 
masses derived from the traditional dynamical method and the gravitational 
lensing based on the published data in literature (Wu 1994; Wu & Fang 1996, 
1997; Wu et al. 1998a). In this talk, I summarize the most recent finding of 
such series. The purpose is to demonstrate how large the uncertainties could be 
among the present various cluster mass estimates. Such a statistical comparison 
may eventually help towards confirming or resolving the baryon crisis and the 
flo disprepancy, if a flat universe is preferred by the standard inflationary model, 
and a better understanding of the dynamical evolution of galaxy clusters. 


2. Sample 

In order to avoid prejudices that may have arisen from our selection criteria, we 
adopt only the published data sets in literature and make no attempt to extrapo¬ 
late the original work. For example, we will not make additional computation of 
the X-ray cluster masses that were not provided by the authors, even if we have 
known the gas density and temperature distributions. Based on this strategy, 
an extensive search for literature yields 21 measurements for galaxy masses in 
clusters, 320 data for mass of intracluster gas, 79 data for virial cluster masses 
derived from optical galaxies, 152+226 data for X-ray cluster masses derived 
from X-ray emitting gas, 8 data for X-ray group masses, and 46/33 data for 
the projected cluster masses derived from strong/weak lensing. Compared to 
the cluster sample of White, Jones and Forman (1997 and hereafter WJF), the 
present sample extends the data of baryon fractions from 176 to 304, the X-ray 
cluster masses from 226 to 378, the gas masses from 226 to 320. Note, however, 
that the data points are apparently larger than the actual cluster population 
because in many cases several measurements have been carried out toward a 
single cluster. We have simply added the observed/derived data that are not 
among the WJF sample. Our data sets and references will be reported elsewhere 
(also available upon request). 

3. Baryon fraction 

Figure |T] displays the 304 measurements of cluster baryon (gas) fractions ft, for 
238 clusters obtained at different radii. With the data being properly binned 
according to radius, an increasing tendency of fj, along cluster radius is clearly 
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Figure 1. (a)The baryon fractions in clusters of galaxies obtained 
from X-ray observations under the assumption of isothermal and hy¬ 
drostatic equilibrium. A total of 304 data points for 238 clusters in 
literature are shown. (b)The same as (a) but the data sets are binned 
such that each bin contains 19 measurements. The observed baryon 
fractions in galaxies and groups of galaxies are also illustrated. The 
dashed lines show the BBN prediction for a flat universe of Qq = 1, 
for which we have already accounted for the uncertainty in the recent 
determination of the primordial Deuterium abundance. 


presented. Such a variation was actually noticed two decades ago and has been 
confirmed by a number of recent observations (e.g. White & Fabian 1995; Ettori 
et al. 1997; David 1997; WJF). Two conclusions can be drawn immediately: 
(l)Because the cluster baryon fraction is representative of the universal value, 
the measurements of the baryon fractions over an ensemble of clusters shown 
in Figure [l] thus imply a cosmological value of = Qb/Q m ~ 30%. This leads 
to < 0.3 combined with the BBN prediction of < 0.08. (2)The increase 
of baryon fraction with radius indicates that a greater fraction of dark matter 
is distributed at small scales than at large scales, in contradiction with the 
conventional point of view. 

4. Statistical comparison of mass estimates 

The original data of the observationally determined cluster masses versus cluster 
radii are shown in Figure ||, in which we have also displayed the projected mean 
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Figure 2. Different mass estimates versus radii for galaxy clusters. 
All the data points are taken from literature. Note that the cluster 
masses derived from gravitational lensing are the projected ones. 


cluster mass required to preduce the quasar-cluster associations detected at the 
smallest cluster radius r = 1.5 Mpc among the four measurements (Wu & Fang 
1996a), and the recent result based on the red galaxy counts behind A1689 
(Taylor et al. 1998). 

The X-ray/optical selected clusters exhibit a large dispersion in their dy¬ 
namical masses due to the difference of richness. On the other hand, the lensing 
clusters are usually the most massive ones at intermediate redshifts, which en¬ 
ables them to act as lenses for distant galaxies. This requires that our comparison 
of different cluster mass estimates can only be made among the hot and massive 
clusters. Examination of the lensing cluster sample (see also Wu &; Fang 1997) 
shows that the mean X-ray temperature and galaxy velocity dispersion of lens¬ 
ing clusters are approximately T = 7.5 keV and a = 1200 km s -1 , respectively. 
We slightly relax these limits and use the criteria of T > 7 keV or a > 1100 
km s -1 to select the X-ray/optical clusters for comparison. This leaves us 49/56 
measurements from the WJF/other samples and 18 data for the virial cluster 
masses M v i r . Furthermore, in order to facilitate a statistical comparison among 
different cluster mass estimates, the selected data sets are properly binned and 
shown in Figure [| We fit the X-ray cluster mass M xray using an isothermal f5 
model with core radius r c for the distributions of hot X-ray emitting gas based 
on hydrostatic equilibrium hypothesis (Figure |j). Our best-fit of the X-ray data 
(WJF + others) for (3 = 2/3 to M xray = 2 (kT / G fj,m v )r 3 / (r 2 +r 2 ) gives r c = 0.25 
Mpc and T = 9.5 keV. Apparently, no significant difference between the mean 
values of M xray and M v i r is detected. Since the gravitational lensing always 
provides a projected gravitating cluster mass mi ens , we need to transform the 
best-fitted mass profiles M xray into the corresponding 2-D masses m xray for the 
purpose of comparison. Alternatively, it should be noted that the current weak 
lensing technique only sets a lower bound on mi ens (r ) within cluster radius r 
(e.g. Fahlman et al. 1994): mi ens (r) > Trr 2 T, crit J^ max (yr)(l - r 2 j r 2 max Y x d\nr, 
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Figure 3. (a)Statistical comparison of cluster masses inside radii r 
derived from dynamical analysis and gravitational lensing. We only 
select the massive X-ray/optical clusters in Figure ||| that satisfy T >7 
keV or a > 1000 km s -1 . The dynamical masses for these selected clus¬ 
ters are properly binned for illustration and fitting. The dashed/solid 
lines are the best-fit of the 3-D/2-D X-ray cluster masses to the theo¬ 
retically expected results under the hydrostatic equilibrium hypothesis 
for an isothermal /3 gas distribution with r c = 0.25 and T = 9.5. Multi¬ 
plying the best-fitted 2-D X-ray cluster mass (solid line) by a factor of 2 
gives the result shown by dotted line. The binned gas mass distribution 
of all the measurements (320) is also shown. 


where X cr j t = (c 2 /4 ttG)(D s /D^D^ s) is the critical mass density, with D^, D s 
and Dfi s being the angular diameter distances to the cluster, to the background 
galaxies, and from the cluster to the galaxies, respectively, ( 7 t ) denotes the 
shear effect on the image configuration of background galaxies introduced by 
the intervening cluster gravitational potential, and r max is the maximum radius 
of a control annulus. 

It turns out that the dynamical masses of clusters given by the traditional 
methods, namely X-ray analysis and virial theorem, are in good agreement with 
the gravitating cluster masses revealed by weak lensing. While we cannot exclude 
the possibility that dynamical analysis may systematically underestimate the 
true cluster masses as compared with strong lensing, the mass discrepancy, if 
any, is well within a factor of ~ 2. When the measurements errors are included 
and the uncertainties in our selection criteria of T > 7 keV or a > 1000 km s _1 
are taken into account, we conclude that there is no significant mass discrepancy 
among the various cluster mass estimators. This is compatible with the previous 
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Figure 4. The a-T relationship for the 141 clusters in our sample. 
The low-redshift (2,0.1) and high-redshift (z > 0.1) clusters are repre¬ 
sented by the open triangles (106) and the filled squares (35), respec¬ 
tively. The solid line is the (unweighted) best-fit to the data. 


similar work made for individual lensing clusters (e.g. Allen et al. 1996; Squires 
et al. 1996; Markevitch 1997; Small et al. 1997; Allen 1997; Wu & Fang 
1996a,b,1997). 

5. The a-T relationship 

The consistency between the different mass estimates indicates that as a whole, 
galaxy clusters can be regarded as the well relaxed dynamical systems, though 
substructures and complex temperature patterns have been observed at small 
scales. A simple and robust way to test this scenario is to employ the velocity 
dispersion - temperature relationship: If both the galaxies and gas are the tracers 
of the depth and shape of a common gravitational potential, we would expect 
a ~ T 0 ' 5 (Cavaliere & Fusco-Femiano 1976). In Figure || we show the updated 
a-T relationship based on 141 clusters selected from our largest cluster sample, 
for which both a and T are observationally determined, where we exclude those 
clusters whose a and/or T are obtained by indirect methods such as the L x - 
a and L x - T correlations. Our best-fitted result reads a = io 2 - 5 ' ±0 - 03 T a53±0 ' 04 
(unweighted) and a = io 2 - 53±0 - 02 T°’ 60±0 ’ 04 (weighted), which are fully consistent 
with previous similar work but with a reduced s.d. of the residuals (see Wu et 
al. 1998b and references therein). Therefore, at 3<r confidence level the a-T 
relationship also supports for cluster of galaxies being a virialized system. 

6. Discussion and conclusions 

The overall observed baryon fraction of galaxy clusters is ~ 30%. This value has 
been justified by our statistical comparison of various mass estimates of clusters 
published so far in literature, which shows a consistency between the cluster 
masses obtained by traditional dynamical methods and gravitational lensing, 
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while the latter is independent of dynamical state and matter content. There¬ 
fore, such a cosmological baryon fraction can be used for the determination of 
the mean mass density of the universe. In combination with the BBN predic¬ 
tion an open universe of Q m < 0.3 is thus preferred. It turns out that the 
solution to the so-called baryon crisis is either to modify the standard BBN 
model or inflationary model, or to introduce a nonzero cosmological constant 
into physics/astrophysics. 

The increasing tendency of the observed cluster baryon fraction to ~ 30% 
with scales (Figure |]) has brought about not only the baryon crisis, if we live in 
a flat universe and if we do not intend to accept a nonzero cosmological constant, 
but also a mystery that a greater fraction of dark matter exists at small scales 
than at large scales. Future cosmological study of structure formation should 
take this fact into account. Meanwhile, we should be aware that the fraction of 
baryonic matter in the universe is not small at all: ~ 30% of the matter of the 
universe is actually visible ! 
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